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Abstract 

Noble metals adopt close-packed structures at ambient pressure and rarely undergo structural transforma- 
tion at high pressures. Platinum (Pt), in particular, is normally considered to be unreactive and is therefore 
not expected to form hydrides under pressure. We predict that platinum hydride (PtH) has a lower enthalpy 
than its constituents solid Pt and molecular hydrogen at pressures above 21.5 GPa. We have calculated 
structural phase transitions from tetragonal to hexagonal close-packed or face-centered cubic (fee) PtH be- 
tween 70 and 80 GPa. Linear response calculations indicate that PtH is a superconductor at these pressures 
with a critical temperature of about 10-25 K. These findings help to shed light on recent observations of 
pressure-induced metallization and superconductivity in hydrogen-rich materials. We show that formation 
of fee metal hydrides under pressure is common among noble metal hydrides and examine the possibility 
of superconductivity in these materials. 
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Pt maintains the face-centered cubic crystal structure up to at least 304 GPa [1], which covers 
most of the range currently accessible in static diamond anvil cell (DAC) experiments, and makes 
it suitable for use as a pressure standard Q2] |3J. Pt is frequently used to form electrodes (which 
allow measurements of electrical resistance, the concentration of electrical carriers and the Hall 
mobility) and as a heat coupler. As a noble metal, Pt is fairly unreactive at ambient conditions, 
and this has generally been assumed also to hold at high pressures. However, chemical reactions 
of Pt with other materials present in the sample chamber can lead to the unexpected formation of 
new compounds as, for example, described by Gregoryanz et al. y4] who observed formation of 
Pt nitride when heating nitrogen with a piece of platinum. 

The hydrogen-dense material silane (SiH 4 ) forms molecular crystals at low pressures. Two ex- 



perimental studies of the metallization of SiH 4 under pressure have recently appeared 051 |6J. The 
first study [|6|] provided optical evidence for metallization of SiH 4 above 60 GPa. Remarkably, the 
second study by Eremets et al. [5] reported superconductivity of SiH 4 in the pressure range 60-150 
GPa with a superconducting transition temperature as high as 17 K. These experimental findings 
have motivated a number of theoretical investigations of compressed SiH 4 [7-9] suggesting that 
metallization and/or superconductivity can be achieved under pressure, but a clear explanation 
of the observations has not yet been achieved. While the above two experimental studies Q5L |6Q 
reported metallization of SiH 4 above 60 GPa, whereas some other experiments claimed decompo- 
sition, amorphization, and the formation other metal-hydrogen alloys besides silane [fiol . 11]. This 
complexity itself is associated with hydrogen as being extremely light and mobile, as well as ex- 
traordinarily reactive. The properties of hydrogen are becoming drastically enhanced under high 
pressure and temperature, and make the interpretation of experimental data for hydrides ||12l-ll4|| 
rather inconclusive. 

Degtyareva et al. Ill ill have speculated that the formation of PtH under high pressures might 
explain some of the observations reported in Ref. 3, and they provided indirect evidence to support 
the hypothesis of PtH formation. Moreover, it appears likely that the silane in the experiments 
of Eremets et al. |5Q partially decomposed, leading to release of hydrogen. Indeed SiH 4 is a 
metastable compound at low pressures which is thermodynamically unstable to decomposition 
into H 2 and Si 111 ill , so that release of hydrogen might occur when the compression is sufficient 
to break bonds in the SiH 4 molecules. Released hydrogen may then have reacted with the Pt 
electrodes to form PtH as the pressure was further increased. 

Confirming the formation of superconducting PtH at high pressures would change the inter- 
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pretation of the data of Eremets et al. [50, but would in no way diminish the importance of the 
experiments. In the new interpretation these experiments will have led to the discoveries of super- 
conducting PtH under pressure and a new insulating phase of SiH 4 , and to the understanding that 
compressed SiH 4 is prone to decomposition. Indeed, Degtyareva et al. [11] have suggested that 
hydrogen release from SiH 4 at high pressures could be useful in synthesising metal hydrides. 

The prevalence of DAC measurements and the steadily growing interest in high-pressure ex- 
perimental studies of metal hydrides raises the possibility of decomposition and hydrogen release 
in such experiments. Inadvertent and undetected formation of noble metal hydrides could lead to 
erroneous pressure estimates and other effects which could substantially alter the interpretation of 
the data. A careful investigation of the possibility of pressure-induced reactions between hydrogen 
and noble metals such as Pt is therefore very important and timely. 

We have used ab initio random structure searching (AIRSS) to predict crystal structures of PtH 
under pressure 111 511 . This approach has been used to predict the structures of materials such as 
hydrogen [16], metal hydrides |n|. transition metals [18 1 , and molecular crystals |19h21|. In 
this study we have used the simplest version of the AIRSS technique in which an ensemble of 
structures is prepared by generating unit cells of random shapes with reasonable volumes into 
which the atoms are placed at random positions. The cell shapes and atomic positions are then 
relaxed to an enthalpy minimum at fixed pressure. We studied simulation cells containing 1, 2, 4, 
6, and 8 formula units of PtH. We found three particularly low -enthalpy metallic structures of PtH 
(the corresponding space groups are given in parentheses): fee (Fm3m), hep (PG-s/mmc), and 
tetragonal (J4m2). 

We have assumed that PtH forms in the 1:1 stoichiometry. The fee structure was used for Pt 
and an hep molecular form for hydrogen [16J, in which the H 2 molecules are randomly oriented in 
space. The variation of the enthalpies H of the fee, hep, and tetragonal phases of PtH relative to 
thatofPt+ ±H 2 is given in Fig.[T] which shows that PtH is more stable at pressures above approxi- 
mately 20 GPa. We also calculated the dynamical stability of the PtH structures and evaluated their 
electron-phonon coupling (EPC) strengths, A, in order to study the possibility of phonon-mediated 
superconductivity in PtH. Since the zero-point motion of the hydrogen nuclei is substantial, we 
have included effects from the zero-point enthalpy (ZPE) calculated within the quasiharmonic ap- 
proximation. For details of the computational methods we refer the reader to the supplementary 
material [22]. It should be noted that there are alternative ways of calculating the superconducting 
transition temperature T c . In particular, there has recently been considerable progress within the 
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framework of density functional theory in developing truly ab initio quantitative predictions of T c 
in the absence of adjustable parameters, which combine electron-phonon coupling and Coulomb 
interactions [23-28]. 
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FIG. 1: (color online). Comparison of the enthalpies of Pt + 5H2 and the tetragonal and fee phases of PtH. 



FigureQ]shows that it is thermodynamically favourable for tetragonal PtH to form from elemen- 
tal Pt and H 2 above 18.7 GPa when the zero-point energy of the nuclei is neglected (above 21.5 
GPa when the zero-point energy contribution is included; see the supplementary material I122ll ). 
Experimentally, a tetragonal PtH was observed at room temperature, 26 GPa on compression and 
20 GPa on decompression Il29ll . We find a transition from tetragonal to hep or fee PtH between 70 
and 80 GPa (Fig.Q]). The fee and hep structures differ only in the stacking of layers, making them 
virtually degenerate in energy (to within a few meV per formula unit). Under pressure, a previous 
calculation using a simple model suggested that the site occupied by hydrogen in <i-metal mono- 
hydrides tends to change from the tetrahedral to the octahedral site and that the hydride transforms 
into the hep or fee structure [30J, in good agreement with our findings. In the present context it 
is more likely for fee PtH to form, since the Pt atoms as the heavier host material are arranged in 
that structure prior to exposure to hydrogen. Experimentally hep PtH not fee PtH was observed 
and we therefore speculate that kinetics might favor a sequence of structural transformation from 
tetragonal PtH to hep PtH. 

Having established the likely formation of PtH from Pt and H 2 in the pressure range relevant 
for the experimental observations reported in Ref. bl we evaluated the electron-phonon coupling 



(EPC) strength in PtH to examine its potential for superconductivity, in a similar fashion to our 
studies of metal tri-hydrides [31, 34], as summarized in the supplementary material I22I . This 
study is crucial because the formation of superconducting SiH 4 was a central claim of Ref. 5, and 
perhaps superconducting PtH could provide an alternative explanation. 
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FIG. 2: (color online). Superconducting transition temperature T c of PtH as a function of pressure. Filled 
squares and hollow hexagons show values of T c for the fee and hep structures, respectively, with the higher 
T c value at each pressure corresponding to the choice /x*=0. 10 and the lower T c value to the choice /i*=0. 13. 
The inset displays the evolution with pressure of the electron-phonon coupling parameter, A, and the asymp- 
totic phonon momentum Qi og (see the supplementary material [22]). The fee phase is dynamically unstable 



below 77 GPa and the hep phase is dynamically unstable below 100 GPa. 



3 apaconstantopoulos suggested that PtH could be superconducting under ambient conditions 



113 511 . but subsequent experiments found no evidence of PtH formation at low pressure Q29 



Tetragonal PtH becomes more stable than Pt + |H 2 above 20 GPa. We calculated the EPC of 
tetragonal PtH in the pressure range 20-60 GPa and found no significant coupling, so it is not 
likely that superconducting tetragonal PtH could exist. Figure [2] shows the variation of the cal- 
culated superconducting transition temperature T c of PtH with pressure. In the whole pressure 
interval of 80-200 GPa we used a fixed value for the effective Coulomb interaction parameter 
H* = 0.13, as has been common procedure [113j,|36|]. However, we also explored the effects on T c 
of changing fi* to 0.10 Ill2l ll3ll. With this, we have covered the two extreme choices of /x*, result- 
ing in corresponding high- and low-estimates of T c (Fig. [2]). It is worth noting that the value of u* 
may also be chosen to vary as a function of the electronic density of states at the Fermi level [39]. 
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The highest value of T c is found in the vicinity of the onset of dynamical stability of the structures, 
as expected for phonon-mediated superconductivity 113 ll — 13311 . The fee structure is dynamically sta- 
ble above 77 GPa where it possesses the highest T c of around 25 K. T c decreases monotonically 
and quite rapidly as the pressure rises, as shown in Fig. [2] and reaches essentially zero around 140 
GPa. The hep structure is dynamically stable above 100 GPa (and is almost degenerate in enthalpy 
with fee PtH) and its maximum T c of around 15 K occurs at 100 GPa. 
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FIG. 3: (color online). Spectral function a 2 F (red shaded areas), the integral of the spectral function up to 
frequency uj ( a 2 F(ui')/uj' dui') (red line), and the phonon density of states (solid black line) as a function 
of frequency at selected pressures for (a) fee PtH and (b) hep PtH. 



Figure|3]shows the spectral function a 2 F and the phonon density of states (PhDOS) at different 
pressures for both fee and hep PtH. At the lowest pressure at which the fee phase is stable (77 GPa), 
the Pt and H atoms both contribute to the EPC (Fig. [3k). With increasing pressure, the spectral 
function splits into high and low frequency bands as shown in the second and third panels of Fig. 
[3^, and only hydrogen vibrations contribute to the EPC in the higher pressure band. The situation 
is similar in the hep phase (Fig. [3})), in the sense that the respective contributions to the EPC 
parameter A from Pt and H are nearly the same as at lower pressures. With increasing pressure, 
the integrated A decreases as the phonon frequencies increase. This behavior can be understood 
from Eq. 3 of the supplementary material D22J, where it can be seen that the spectral function a 2 F 
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is divided by the phonon frequencies. 

Interestingly, Eremets et al. reported an apparent divergence in T c between 80 and 100 GPa |5|] 
which coincides with our calculated transition pressure between fee and hep PtH. These structural 
transformations are accompanied by softening of the acoustic branch modes which is a precursor 
of structural transformation. High-pressure experiments on PtH have been reported earlier [29J, 
but since they were performed at room temperature, no information regarding possible supercon- 
ductivity could be obtained. It would be very interesting to study compressed PtH at temperatures 
below 25 K to determine whether it is indeed superconducting. 

At low pressures, fcc-Pt absorbs hydrogen and forms tetragonal PtH, which is a distortion of the 
fee structure. Underpressure, the tetragonal PtH phase is restored to its original fee structure or the 
similar close-packed hep structure, which possess strong EPC and hence a high superconducting 
critical temperature. Such a structural transformation from a lower-symmetry phase to the fee 
structure also applies to other hydrides formed by metals which are neighbors of Pt in the periodic 
table (referred to in the following as the Pt family of hydrides). We investigated the metals Rh, 
Pd, Ag, Ir, and Au, which all adopt the fee structure at ambient pressure and have extremely small 
or zero superconducting T c , and determined the onset pressure at which the corresponding fee 
hydrides become dynamically stable. We find that fcc-RhH (lattice constant a = 4.02 A) and fec- 
PdH(a = 4.08 A) are stable at almost ambient pressure, consistent with experimental observations 



42, 



43 



4511 . Also, our calculations are in good agreement with the measured lattice constants 



which were reported to be 4.01 A at 3.8 GPa (hydrogen gas pressure) for fcc-RhH Q42 



43D and 



4.09 A at ambient pressure for PdH D45|]. The metal hydrides fcc-AgH (a = 4.00 A) and fcc-IrH (a 
= 3.88 A) are dynamically stable at 50 and 80 GPa, respectively. Dynamical stability of fcc-AuH 
(a = 3.80 A) is achieved only above 220 GPa. We emphasize that other hydride structures besides 
fee might be stable at lower pressures. 

Similarly to PtH, the dynamical stability of other fec-structure hydrides can be regarded as 
precursors for superconductivity, while lower symmetric structures might form at much lower 
pressures [|44[|. Figure |4] shows the calculated T c at the onset of dynamical stability of the fee 
structure in the Pt family of hydrides. PdH is calculated to possess a T c of 20 K, which is com- 
patible with experiments on non-stoichiometric PdH a . which reported a T c of 9 K 14511 . although 
we note that our results do not reproduce the observed isotope dependence of T c in PdH/PdD [14]. 
We predict IrH to have a T c of 7 K at 80 GPa. Remarkably, the T c of fcc-AuH is estimated to be 
21 K, however, it takes a very high pressure of 220 GPa to reach this superconducting state. As 
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FIG. 4: (color online). Maximum T c (calculated at the onset of dynamical stability of the fee phase) for the 
metal hydrides RhH, PdH, AgH, IrH, PtH, and AuH. Note that the corresponding elemental metals in their 
fee form all possess either extremely small T c (Rh and Ir) or are not superconducting at all (Pd, Ag, Pt, and 
Au). 

the pressure is increased, the T c of all these fec-hydrides will decrease, analogously to the case of 
PtH. We predict that AgH is not a superconductor, because the electronic density of states exhibits 
a small bandgap at the Fermi level. AgH does however possesses the highest Cli og , and therefore, 
if electrons were to be donated in a complex alloy, this could lead to a good superconducting ma- 
terial, such as (PdAg)H;,, possessing a T c of about 16 K [47]. RhH has a small T c and it is worth 
note that no superconductivity of RhH has been detected at T c > 0.3 K 1480 . 

In conclusion, we have shown from first principles that formation of PtH from Pt and H 2 can 
occur underpressure. Based on the empirical Allen-Dynes equation (see supplementary material), 
the resulting product is predicted to be a superconductor. Throughout the studied pressure interval 
of 80-200 GPa, we used the effective Coulomb interaction parameter ji* — 0.13, a value which has 



often been used [36]. We did, however, consider the effect of changing /i* to 0.10 on the predicted 
T c . If our prediction concerning the formation of superconducting PtH at high pressures would be 
experimentally confirmed then this might lead to new interpretation of the data of Eremets et al. 
|5|]. That work would mark the discovery of superconducting PtH under pressure and a new insu- 
lating phase of SiH 4 , and to the understanding that compressed SiH 4 is prone to decomposition. 
We expect these findings to be relevant to recent observations of pressure-induced metallization 
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and superconductivity in hydrogen-rich materials. Indeed, our calculations show that the forma- 
tion of superconducting fee metal hydrides underpressure may be a common phenomenon among 
noble metal hydrides. The pressure-induced formation of PtH and of other noble metal hydrides, 
along with their associated superconductivity, represents an important aspect which should be con- 
sidered in the interpretation of diamond anvil cell experiments whenever the corresponding metals 
and hydrogen are present, with the latter introduced either directly or entering via decomposition 
of hydrides. 
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Supplementary Material 



"Formation and superconductivity of platinum hydride under pressure" 
D. Y. Kim, R. H. Scheicher, C. J. Pickard, R. J. Needs, and R. Ahuja 



A. Methods 

Our results were obtained using density-functional-theory (DFT) methods and the 

generalized gradient approximation (GGA) of Perdew, Burke and Ernzerhof for the exchange- 
correlation functional Q3L |4£| . The structure searching was carried out with the CASTEP plane- 
wave code [5] and ultrasoft pseudopotentials [g]. For the searches we used a plane- wave basis-set 
cutoff of 550 eV and a Brillouin-zone integration grid of spacing 2%x 0.07 A -1 . The structures of 
interest were further relaxed for enthalpy comparisons, which were done with the VASP code [7, §| 
with the projector-augmented wave (PAW) approach [9]. A cutoff energy of 700 eV was found 
to produce converged results. The Monkhorst-Pack scheme ifioll was used to generate the Appoint 
grid meshes of 31 x 31 x 31 for Pt and PtH (T-centered in the case of hcp-PtH). For calculations of 
solid hydrogen, a large supercell containing 36 H 2 molecules was used; therefore, the Brillouin- 
zone was only sampled at the r -point. Calculations of the phonon modes and frequencies were 
perf ormed with density-functional perturbation theory using the Quantum ESPRESSO code 
i t 1 1|1 - The electronic orbitals were expanded in a plane- wave basis set with a kinetic energy cutoff 
of 60 Ry. The Brillouin zone (BZ) integrations in the electronic and phonon calculations were 
performed using Monkhorst-Pack meshes. We refer to meshes of Appoints for electronic states and 
meshes of g-points for phonons. The EPC matrix elements were computed in the first BZ on a 
8x8x8 g-mesh (8x8x4 g-mesh) using individual EPC matrices obtained with a 24 x 24 x 24 
(24 x 24 x 16) fc-points mesh for fee (hep). Convergence tests of the EPC calculations were 
performed with uptol2xl2xl2 g-meshes for fee and 8x8x8 g-meshes for hep. 

The total phonon zero-point energy was calculated as E p h = \ J2i over a range of vol- 
umes V, where hwi are the individual phonon energies. We fitted the results to a polynomial and 
differentiated with respect to V to obtain the zero-point enthalpy as H ph = E ph — VdE ph /dV . 

T c was calculated using the Allen-Dynes equation. J12I1 

T-^expf 1-04(1 + A) \ 

lc ~ 1.2 P V A-/i*(l + 0.62A);' W 
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with the standard value //* = 0.13. The normalized weighting function of the Eliashberg theory 
llnll is 

<?M = ^ 2 F{u). (2) 
The parameter A is a dimensionless measure of the strength of a 2 F: 

POO 

X = 2 dua 2 F(u)/u, (3) 
Jo 

and the logarithmic average frequency, Vti og in units of K, is 

Q log = exp g(tu) Inu doj^j . (4) 
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B. Crystal structure details 



tetragonal PtH at 46.6 GPa 

lattice vectors: 

a = b = 3.422 A, c = 2.782 A 

atomic positions: 
Pt 1/2 3/4 

Pt 1/2 1/4 

H 

H 1/2 1/2 1/2 

fee PtH at 40.6 GPa 

lattice vector: 
a = 2.843 A 
atomic positions: 
Pt 1/2 1/2 1/2 

H 

hep PtH at 48.1 GPa 

lattice vectors: 

a = 2.789 A, c = 4.733 A 

atomic positions: 
Pt 1/3 2/3 1/4 

Pt 2/3 1/3 3/4 

H 

HO 1/2 



C. Phonon enthalpy contribution 



Figure|5]shows the effect of the nuclear zero-point enthalpy (ZPE) H p h, as defined in the Meth- 
ods section, on the formation pressure of PtH. The difference in total enthalpy between tetragonal 
PtH and Pt + | H 2 is plotted. As can be seen, the inclusion of H ph shifts the formation pressure 
up to 51.8 GPa, far beyond the 18.7 GPa obtained when the ZPE is neglected. Such a shift seems 
unreasonably high, and we therefore decided to calculate the effect of H ph for a material with a 
well-established formation pressure. 

PdH is known to form close to ambient pressure [U4I1 . Figure |6] shows that the enthalpy compar- 
ison of PdH with Pd + | H 2 yields a formation pressure in the vicinity of GPa, actually slightly 
below it. However, once H p h is included, the formation pressure shifts to 23 GPa (Fig. |6]), in 
contradiction with experiment. 

What is the source of this discrepancy? The zero-point enthalpy H ph contributions from the 
metal and metal hydride can be rather precisely determined computationally, in particular because 
there is very little uncertainty about the slope dE ph /dV, the phonon pressure term, which is mul- 
tiplied by the volume and added to the phonon energy to yield the phonon enthalpy H ph . For Pt 
and PtH, the phonon energy E p h is found to depend linearly on the cell volume, and so dE p h/dV 
remains virtually constant throughout the pressure range studied. The main uncertainty stems 
from the contribution to H ph from H 2 , a system which is notoriously difficult to describe properly 
within density functional theory at low pressures 111 511 . As a consequence, it is hard to determine 
the phonon energy E p h accurately at low pressures, and furthermore, the non-linear dependence of 
E ph on V in solid hydrogen (see Figure [7]) prevents us from obtaining a reliable value of dE ph /dV 
at lower pressure. 

To correct for the error in E ph and in the slope dE ph /dV, we therefore determined a correction 
factor which multiplies H ph for H 2 , in order to recover the known formation pressure of PdH. 
Figure |6] shows the enthalpy when H ph for H 2 is scaled by a factor of 3.3, which yields a reasonable 
formation pressure in the vicinity of GPa. 

Applying the same scaling correction to H p h for H 2 in the enthalpy comparison for PtH (Fig. 
[8]), we obtain a formation pressure of 21.5 GPa, corresponding to a much more reasonable shift of 
2.8 GPa due to the effects of ZPE. 
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PdH 

-. Pd + 1/2 H 2 

_ PdH (with H ph ) 

Pd + 1/2 H 2 (with original H ph ) 

Pd + 1/2 H 2 (with corrected H ph ) 

i i i i I i I i I i 

10 20 30 40 50 

Pressure (GPa) 

FIG. 6: Comparison of the enthalpies of Pd + i H2 and PdH without H p h, with H p h, and with the corrected 
Hph contribution. 
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FIG. 7: Eph for solid H2 is plotted as a function of volume. The calculated data are shown as black circles. 
The connecting red lines are merely a guide to the eye. The vertical dashed line indicates the equilibrium 
volume of hydrogen at zero pressure. 
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2 ' 1 ' 1 ' 1 ' 1 ' 1 ' 

'0 20 40 60 80 100 

Pressure (GPa) 



FIG. 8: Comparison of the enthalpies of Pt + i H2 and tetragonal-PtH, including the correction discussed 
in the text. 
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